A comparison of various measures for the average transport coefficient q 
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Jet quenching, i.e. the suppression of high transverse momentum (Pt) hadron production in 
ultrarelativistic heavy-ion collisions is among the most striking experimental signatures of bulk 
medium formation. Efforts with the aim of extracting quantitative information about the bulk 
medium from the measured suppression mainly focus on the extraction of an averaged transport 
coefficient (q) as a measure of the medium jet quenching power, with the underlying assumption 
that (q) is a meaningful quantity to make comparisons both among different models and between 
models and data. In this note, the main uncertainties associated with the extraction of q from 
model fits to data are briefly reviewed before the notion of an average transport coefficient itself is 
investigated. It is shown in a case study that the choice of a meaningful average is far from unique 
and that different well-motivated definitions of averages lead not only to differences in the outcome 
of about a factor four, but also in a different change of (<J) with collision centrality. This casts 
some doubt on the idea that condensing the information about jet-medium interactions into a single 
average parameter is a meaningful procedure. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The suppression of high transverse momentum (Pt) 
hadron production in heavy-ion (A-A) collisions with re- 
spect to the scaled expectation value from p-p collisions is 
experimentally well established from both single hadron 
observables [J, [|| and hard dihadron correlations 0,0. 
Theoretical models using well-constrained fluid dynami- 
cal simulations of the bulk matter evolution as the back- 
ground in which high px processes take place can describe 
these observables qualitatively well 0] , thus the natu- 
ral next step is to extract quantitative parameters of the 
medium which characterize its ability to modify high Pt 
partons from the experiments. 

Since detailed differential observables like the average 
medium induced parton energy loss AE probability 
distribution P(AE) are difficult to obtain from single 
hadron observables as the nuclear suppression factor Raa 
is insensitive to the functional form of P(AE) @,H|, the 
focus in recent times has been on the determination of 
the averaged transport coefficient (q) (the average virtu- 
ally gain per unit pathlength of a parton travelling in the 
medium) or an equivalent parameter within various mod- 
els of energy loss from detailed fits to single hadron and 
dihadron observables (see e.g. H, E3] for such attempts). 
Unfortunately, there is a history of large differences in the 
value of (q) extracted from different models, and some of 
these differences persist even in a comprehensive analysis 
of three differet models under the same set of assumptions 
about the bulk medium evolution [f| . A couple of reasons 
why different calculations may result in different medium 
parameters have been identified since: 

• The assumptions about the medium geometry and 
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evolution may not be the same. While for ex- 
ample early computations in the Armesto-Salgado- 
Wiedemann (ASW) formalism [ll[ were done in a 
static cylindrical geometry and found an averaged 
q = 5.. 15 GeV 2 /fm, later computations utilized the 
same formalism in a fluid-dynamical background 
(if! EH and found averaged transport coefficients 
almost an order of magnitude smaller. As results 
presented in [13| indicate, even the variation be- 
tween dynamical medium evolution models which 
are constrained by the medium induces an uncert- 
inty if ~ 50 % on the extraction of a q. 
The model connecting thermodynamical parame- 
ters with a value of q may not be the same. Com- 
mon prescriptions in the absence of flow use q ~ s, 
q~T 3 or q ~ e ' where s, T and e are the entropy 
density, temperature and energy density as taken 
from a hydrodynamical description of the medium. 
While all these relations lead to the same result if 
the medium can be described as an ideal gas, they 
result in an additional ~ 50 % uncertainty in the 
extraction of q if e.g. a bag model equation of state 
is used for the medium [5j. 

For formalisms derived using a collinear approxima- 
tion, prescriptions to treat finite kinematics must 
be used in the application of the formalism to de- 
scribe data. In doing so, a large systematic uncer- 
tainty is introduced by the detailed choice of the 
finite kinematics prescription [lij . 



II. AVERAGING PROCEDURES 

However, even if all these issues can be overcome, any ex- 
traction of a q from an evolving medium eventually faces 
the question how a meaningful average can be defined. 
From a theoretical perspective, q is a local property of 
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the medium. In the static case, it is tied to local medium 
thermodynamics, leading to the dependence q(r) s ,r, <f>, r) 
on the complete set of parameters spacetime rapidity r) s , 
radius r, angle 4> and proper evolution time r. In a flow- 
ing medium, q acquires a dependence on the relative ori- 
entation of the direction of parton propagation and the 
flow velocity as [l5[ 



Qstatic (cosh p - sinh p cos a) 



(1) 



where p is the flow rapidity and a the relative angle be- 
tween flow and parton trajectory. In principle, the value 
of q should also depend on the resolution scale E at which 
the parton probes the medium (l6| . Trying to condense 
the full information contained in q(ji s , r, </>, r, p, a, E) into 
a single number is a bit like asking a Meteorologist to 
summarize the weather forecast for a place for the next 
48 hours in a single number — it may not nesessarily 
be a meaningful exercise. In the formalism comparison 
study the problem has been circumvented by com- 
paring every energy loss model for the same hydrody- 
namical medium and making the single highest value of 
q the standard of comparison. This works as almost the 
whole parameter dependence of q is specified by the hy- 
drodynamical medium, but results in values of q which 
are only meaningful for this particular evolution model 
and do now allow a meaningful comparison to results ob- 
tained in other models for the bulk medium. 
In the following, several different averaging procedures 
are discussed. Note that all computations are done for 
the same ASW energy loss formalism and the same 3- 
d hydrodynamical bulk matter evolution model as de- 
scribed in detail in In all cases, the model repro- 
duces Raa(Pt) for centeral and non-central collisions. 
Any differences in the extracted values of (q) quoted be- 
low are thus exclusively caused by the different averaging 
procedures. 



A. Spacetime averaging 

One possibility to define an average transport coefficient 
is to make averages over the spacetime evolution of the 
hydrodynamical medium, given a model for the relation 
between q and thermodynamical parameters. Dropping 
the flow correction Eq. {!]), a relation like 



q(r) s ,r,<f>,T) = K ■ e(r) s , r, (/>, r) 3/4 



(2) 



can be used to compute the nuclear suppression factor 
Raa as outlined e.g. in 0, 0, with K adjusted such 
that the data is described. For the average transport 
coefficient, the expression 



<<?) = 



/ dr] s drd4>dTq(ri s ,r, <j>,r) 
J di] s drd(j)dT 



where the integration boundaries are given by the ther- 
malization and freeze-out hypersurfaces of the medium. 
If this is done for central collisions at 200 AGeV using 
the hydro model of Q3, a value of (q) = 0.82 GeV 2 /fm 
is obtained. Note that this result is based on exactly the 
same ASW formalism described in Q with a q ma x = 18.5 
GeV 2 /fm. 

However, this prescription has an unwanted dependence 
on the choice of the freeze-out conditions. Increasing the 
freeze-out temperature T F from 130 MeV (the standard 
choice in [l7|) to 140 MeV does not alter the value of K 
obtained in a fit to Raa, but the value of (q) changes to 
1.02 GeV 2 /fm. The physical reason for this is of course 
that hard partons tend to escape from the medium long 
before bulk freeze-out occurs, thus whatever happens at 
late times in the evolution does not influence hard parton 
production, yet the late time behaviour of the hydrody- 
namics still enters the averaging procedure Eq. ([3]). In 
particular, this implies that two models in which the in- 
teraction between hard partons and the medium happens 
in exactly the same way may still lead to different values 
of (q) if their late-time dynamics is different. 
If one cuts the averaging after a finite time with the ar- 
gument that partons only see the evolution before this 
cut, the extracted average transport coefficient is signif- 
icantly larger. Assuming all partons escape by the time 
of 5 fm/c increases (q) = 1.8 GeV 2 /fm, changing the cut 
to 3 fm/c yields 2.33 GeV 2 /fm. While such a cut re- 
moves the artificial smallness of the extracted transport 
coefficient due to the late time evolution of the medium, 
(q) is then dependent on the relatively arbitrary choice 
of a timescale. Thus, spacetime averages do not provide 
a satisfactory solution to all aspects of the problem. 



B. Path averaging 

Since the weakness of the spacetime averaging approach 
is that the fact that hard partons leave the medium after 
a finite time is not taken into account, one may think of 
averaging along a single parton path as 



(q) P a 



th 



(4) 



(3) 



where £ is a variable describing the spacetime trajectory 
of the parton through the medium and the integration 
boundaries are again given by the intersection of the par- 
ton trajectory with the thermalization and freeze-out hy- 
persurfaces. If one assumes that the initial vertices of the 
hard processes responsible for parton production scale as 
the number of binary collision with the overlap 



P / v T A (vo + h/2)T A (v Q ~h/2) 
P{x ,ya) = Tf, — tu , (5) 

where the thickness function is given in terms of 
Woods-Saxon the nuclear density Pa(y, z) as Ta{v) — 
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J dzpA(r,z), the average transport coefficient can be 
computed by averaging Eq. Q over all possible initial po- 
sitions (xq, uq) and angles with a weight given by Eq. ([5]). 
If this is done in the model outlined above, one finds 
(q) = 2.2 GeV 2 /fm. This is consistent with the space- 
time average with a cut in r between 3 and 4 fm/c and 
indicates that the bulk of hard partons leaves the medium 
within this timescale. As expected, the value obtained in 
the path-based averaging is significantly higher than the 
spacetime average over the whole medium, as the path 
averaging shifts weight to the dense medium center and 
to early times. 

However, once one tracks individual parton trajectories, 
one may pose the question if one should take into account 
the fact that propagating gluons carry a different color 
charge and couple by a factor 9/4 more strongly to the 
medium. If the relative weight of quarks and gluons is ob- 
tained in a leanding order perturbative Quantum Chro- 
modynamics Calculation, a value of (q) — 3.4 GeV 2 /fm 
is obtained with a small Pt dependence caused by the 
change in relative quark and gluon yield at different Pt- 
But one may also argue that averaging over all parton 
trajectories is not suitable to make a comparison with 
the data, because the measurement of the nuclear sup- 
pression factor reflects only partons which, after energy 
loss and fragmentation, have produced a sufficiently high 
Pt hadron. If one supplements the path-based averaging 
with a condition that the leading hadron after fragmen- 
tation must be found above a certain Pt scale, the av- 
erage weights change somewhat and the result acquires 
a scale dependence (which is in practice small for RHIC 
kinematics but may be noticeable at the LHC). For a 
scale of 6 GeV, (q) — 2.1 GeV 2 /fm is found. This num- 
ber is smaller than the one found when averaging using 
all paths, as the trigger condition removes weight from 
the dense regions of the medium where partons are ab- 
sorbed. In extreme cases, this may be a problem: If one 
assumes a threshold core-halo scenario in which partons 
are either absorbed by a dense core or escape undisturbed 
from a halo, the path-based averaging with trigger condi- 
tion would find (q) = since only partons from the halo 
without medium effect are observed. 



C. Centrality dependence 

Conceptually, one wants to define (q) as a quantity which 
is extracted from and related to a measured quantity such 
as Raa-i but on the other hand characterizes the pro- 
duced medium. Therefore, any (q) should reflect a change 
of the medium properties in a characteristic way. This 
can for example be tested by studying the relative change 
of (q) when going to non-central collisions in the hydro- 
dynamical evolution. Fig. [T] shows the average transport 
coefficient, normalized by its value for b — 2.4 fm, for the 
various averaging procedures described above along with 
the maximal q measure as used in [j| when the medium 
is changed to non-central collisions. 
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FIG. 1: (Color online) 



It is obvious from the figure that the relative change of (q) 
with centrality docs not agree among different averaging 
procedures. The underlying reason is again that differ- 
ent regions within the medium volume receive different 
weight factors in different averaging procedures. 



III. DISCUSSION 



One can roughly classify the various measures of (q) 
according to how much they focus on medium proper- 
ties as opposed to the measurement process in which 
the medium is probed. Spacetime averages of the hy- 
dro medium incorporate the full information about the 
medium, but none about the physics underlying Raa- 
Consequently, they provide a good standard of compar- 
ison among models of the bulk medium, but given the 
constraint that models should describe the experimental 
data for Raa, different bulk models may lead to large 
variations in the extraction of (q) as the averaging probe 
spacetime volumes of the medium which are never probed 
by hard processes. 

In contrast, averages based on actual parton propagation 
and energy loss tie different bulk evolution models much 
closer to the measurement. But for these averaging pro- 
cedures, there are dense regions to which the averaging 
process is blind as almost no partons escape. This im- 
plies that two models required to describe the measured 
Raa may result in a similar (q) in spite of the fact that 
the models, e.g. in terms of central energy densities, are 
very different. Thus, constraining a path-averaged (q) 
does not necessarily constrain the underlying bulk evolu- 
tion model well. 

Given these findings, it would appear that a more dif- 
ferential approach to characterize medium properties is 
needed to do quantitative measurements of medium prop- 
erties. 
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